The Distinguishability of Interacting Dark Energy from Modified Gravity 
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We study the observational viability of coupled quintessence models with their expansion and 
growth histories matched to modified gravity cosmologies. We find that for a Dvali-Gabadadze- 
Porrati model which has been fitted to observations, the matched interacting dark energy models 
are observationally disfavoured. We also study the distinguishability of interacting dark energy 
models matched to scalar-tensor theory cosmologies and show that it is not always possible to find 
a physical interacting dark energy model which shares their expansion and growth histories. 



I. INTRODUCTION 

The Universe appears to be undergoing a late-time ac- 
celerated expansion P, @|. A model which includes a 
cosmological constant A and cold dark matter (CDM) 
evolving according to Einstein's theory of General Rel- 
ativity (GR) provides the best description of this Q. 
There are many alternative explanations however, the 
two main classes of which are modified gravity (MG), 
(see 0-01 for reviews), and dark energy (DE), (see 0-@| 
for reviews), and we must rely on observations to dis- 
criminate between them [Icj . 

It is always possible to find a DE model with a time 
varying equation of state parameter which produces a 
given expansion history [Tll - [l3| . so in a worst-case sce- 
nario a DE model could exactly mimick a MG model's 
expansion history, making them indistinguishable. To 
break this degeneracy it is necessary to take differences 
in the growth of structure into account and a great deal 
of effort has gone into distinguishing DE from MG [14 - 
|42| . It has been argued that by finetuning the properties 
of a DE model its structure growth can also be made to 
mimick that of a given MG theory (43l - l45l |. but by em- 
ploying suitable combinations of observables consistency 
tests can be made which should be able to distinguish 
between realistic models [46l |47|. 

The above works focus on minimally coupled DE but 
it's also possible to match the growth and expansion his- 
tories of MG with interacting dark energy (IDE) mod- 
els [48j], (for recent IDE works see [49| and references 
therein). IDE models can look like modifications of 
GR dO], but they should deviate from GR+ACDM in 
a way which is distinct to that of MG [5l|. In this pa- 
per we investigate their distinguishability by testing the 
observational viability of IDE models with their growth 
and expansion histories matched to MG cosmologies, re- 
stricting ourselves to a flat spacetime in the Newtonian 
regime. 

Section [TTJ revisits an example Dvali-Gabadadze- 
Porrati (DGP) [52] model from [48j to examine the obser- 
vational distinguishability of matched IDE/DGP models. 
Section Mil extends the matching procedure used for the 
DGP case to a more general scalar-tensor theory (STT) 
model and again considers whether the matched IDE 
models can be distinguished from their MG counterparts 



observationally. Our conclusions are then drawn in sec- 
tion nvi 
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INTERACTING DARK ENERGY MATCHED 
TO A DGP COSMOLOGY 



For a scalar field model of IDE a general action may 
be written as, 
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where S rn (g fJiV , ip, tp) is the matter action, with ip being 
the matter field. In (4|| the authors matched a gener- 
alised IDE model to a particular choice of DGP model 
which had been fitted to observations. They used the 
IDE potential and coupling functions to match the DGP 
expansion and growth histories respectively. Essentially 
the evolution of the background CDM density p in the 
IDE model is determined by the matching of its pertur- 
bation 8 = 8p/p to that of the DGP model. 

Fig. [1] shows the evolution of 8 and the different back- 
ground evolutions of the IDE and DGP density param- 
eters 51 and ft, (tildes denote MG quantities through- 
out), where Q = 8nGa 2 p / (3H 2 ) with U = a^da/dr 
and r is conformal time. Also plotted for comparison is 
a GR+ACDM model chosen to give O ~ 0.227, (sub- 
script 0's denote present day quantities throughout), in 
line with recent constraints [3j. 

The original example had initial conditions set early 
in the matter dominated era with an initial DGP energy 
density parameter fij « 1, (subscript i's denote initial 
values). The initial IDE density parameter was Qi = 
0.995 and in addition to this solution we plot the result 
of choosing Q { = 0.996 and fl t = 0.997 in Fig. Q] but 
find that there are no solutions with £li >= 0.997, (see 
Appendix Ia]) . 

This means that there is a limit on how closely one 
can hope to match the evolution of the IDE/DGP densi- 
ties through the choice of the boundary conditions on f2. 
This difference should be evident in any quantity which 
depends on the CDM density, for example the sum of the 
metric potentials. 
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FIG. 1: Evolution of the density perturbation (left) and the density parameters (right) for the matched DGP/IDE models, 
each with a different and a GR+ACDM model. 



The perturbed metric in the Newtonian regime may be 
written, 



ds 2 



a 2 [-(l + 2^)dr 2 + (1 - 2<^)j ij dx i dx j ] 



(2) 



Both the DGP and IDE models obey the same evolution 
equation for the sum of the metric potentials ^ and $, 

<9 J <9 4 (1'+$) = 8TrGa 2 p6, d l d l (t> + $) = 8TrGa 2 p5. (3) 

This quantity is plotted in the left-hand panel of Fig. [2] 
as a function of redshift at late times, making clear the 
significant distinction arising between the IDE and DGP 
models from the restriction on the boundary conditions 
for n. 

One way to test for this difference observationally is to 
use the Eg prameter [54j defined by, 



Er 



d z d l (y + <f>) 



(4) 



where in the Newtonian regime 9 = —8', (primes denote 
derivatives with respect to N = ln(a) throughout). Note 
however that this relation does not hold for all IDE mod- 
els, eg. [49(. The numerator in Eq. ([4]) can be measured 
from weak lensing observations, while the denominator 
can be found from peculiar velocity measurements and 
for the models studied here we have, 
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The right-hand panel of Fig. [5] shows Eq at late times 
for the DGP, IDE and GR+ACDM models along with 
some recent observational constraints [53[. We can see 
that although the DGP model is a good fit, even the 
worst-case IDE model with boundary conditions as close 
as possible to those of the DGP model is disfavoured by 
observations. 



III. INTERACTING DARK ENERGY 
MATCHED TO SCALAR-TENSOR THEORY 
COSMOLOGIES 



In the same vein as the previous section, we now ap- 
ply a similar method to a simple STT model in order to 
explore the potential distinguishability for particular pa- 
rameter values. We assume it to be the large scale limit 
of MG models to which local constraints on the gravity 
theory [55| do not apply due to a screening mechanism 
such as the chameleon |56l458| , thus allowing the effect of 
baryons to be neglected. 

The action for a STT model may be written as, 



S, 
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d 4 x. 



16ttG, 



<j>R- -(V0) 2 - U 



(7) 



In general uj and U can be functions of 4> but for our 
purposes we take them to be constant. The acceleration, 
scalar field and density perturbation equations derived 



3 



y + o 
^ + o 





GR+ACDM 

Q = 0.995 

Q 




FIG. 2: Evolution of the sum of the metric potentials normalised at the present day (left) and the Eg parameter (right) 
for the matched DGP/IDE models, each with a different Hi, and a GR+ACDM model. The observational measurement is 
E G = 0.39 ± 0.06(1<t) at an effective redshift of z = 0.3 HI]. 



from this action are. 



Eq. (|M|) now leads to a differential equation for f2, 
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These equations determine the expansion and growth his- 
tories for both the STT and matched IDE models. The 
action for the IDE model, Eq. ([T]), leads to its fluid, scalar 
field, Friedmann and density perturbation equations, 
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5" = - 1 
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where C is the DE/CDM coupling function and V is the 
scalar field potential, both of which are taken to be free 
functions. Using Eq's (|1HI13I) and comparing Eq. (fTUl) to 
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where, 



C = 16ttG 



1 -L QL -L 9^- 
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Eq. (fT5|) is quadratic in fi' and so we choose the root 
which is typically negative initially, (the alternative 
branch typically leads to increasing D, and the limits de- 
scribed in the Appendices are reached before the present 
day). We can now solve Eq's ([5 HTU)) numerically, along 
with the root of Eq. l[Tg]). to find Q, H, S, 5', (f> and <j)' 
at any given N. In this way the freedom in the coupling 
function G is explicitly used to match the evolutions of 
the <5's, while the freedom in the scalar field potential V 
is used implicitly to match the expansion histories via 
the IDE Friedmann constraint, Eq. ([T3")) . The initial con- 
ditions used are, 



Ni -- 
Hi 



-7, 
1, 



1. 



0. 



5' 



(17) 



Qi is chosen so that f2o is the same as the previously 
mentioned GR+ACDM model's present day density pa- 
rameter when oj —¥ oo, while Ui is determined by the 
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FIG. 3: Evolution of the density perturbation (left) and the density parameters (right) for the matched STT/IDE models and 
a GR+ACDM model. The models STT11, STT20 and STT100 have u = 11, u = 20 and u = 100 respectively, with IDE11, 
IDE20 and IDE100 being their matched IDE counterparts. Note that including <j> in the definition of f2 would bring the STT 
models on the right much closer to the GR+ACDM Q curve. Even then however fi > 1 is still possible despite the spacetime 
being flat because the sum of the gravitational scalar field terms in the STT 'Friedmann' equation can be negative. 



choice of Oj due to the 'Friedmann' constraint, 



Ui = 



6Hi 



-fii 



(18) 



As in the case of the earlier DGP example there is a 
limit on how close f2 j can be to fij (see Appendix [Bj . 
Fig. [3] shows results for three different values of uj where 
in each case f2j has been chosen to be as close as possible 
to fij in the spirit of representing a worst-case scenario 
for distinguishing between the IDE/STT models. The 
evolution equation for the sum of the metric potentials 
in the STT model is, 



did* + $) 



where G = Gcf>o, leading to, 



8uG 



n s 



a pS 



(19) 



(20) 



with the IDE expression as before in Eq. ([6]). Fig. 0] 
shows \t + $ and Eq as functions of z at late times for 
the STT models and their matched IDE counterparts. 
Once again the IDE models lie much farther from the 
GR+ACDM case than their MG counterparts, with all 
but that matched to the lo — 100 STT model lying out- 
side of the observational constraints on Eq. 

In [5^| it was shown that constraints on STT models 
from cosmic microwave background, matter power spec- 
trum and local gravity measurements could be avoided 



using a chameleon mechanism, leading to only a weak 
bound of lo > —1.28. Our model here is essentialy a 
Brans-Dicke theory [60j plus a cosmological constant, 
for which lower bounds of lo > 120(2q-) |6ll| and lo > 
97.8 (2cr ) [62j have been found, (Note that [63] give a lower 
bound of lo > 1000(2<r), but see discussions in [6ll |62|). 

The addition of supernova data would significantly im- 
prove constraints on lo, but account would need to be 
taken of local [64| and temporal [65| variation in the 
gravitational scalar field <fi. In [6(| a recovery of GR at 
late-times sufficient to allow the use of supernova data 
was assumed and bounds of lo > 500 — 1000 from future 
data were forecast. If these constraints can be acheived it 
will not be possible to distinguish between our matched 
STT/IDE models with the Eq results we use here, al- 
though with new data of course the Eq constraints could 
also be tightened. 



IV. CONCLUSIONS 

We have shown that although it is possible to construct 
an IDE model which matches the growth and expansion 
histories of a DGP model fitted to observations, even in 
the worst-case scenario, where their density evolutions 
are as close as theoretically possible, the matched IDE 
model can be distinguished by observations. 

For our simple STT model and its matched IDE coun- 
terpart we have calculated a limit on how similar the 
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FIG. 4: Evolution of the sum of the metric potentials normalised at the present day (left) and the Eg parameter (right) for 
the matched STT/IDE models, each with a different ft, and a GR+ACDM model. The models STT11, STT20 and STT100 
have u) = 11, u) = 20 and lu — 100 respectively, with IDE11, IDE20 and IDE100 being their matched IDE counterparts. The 
observational measurement is Eg = 0.39 ± 0.06(lcr) at an effective redshift of z = 0.3 [53 |. 



initial matter densities can be. This limit depends on 
the strength of deviation from GR and we find that in 
cases which differ significantly from GR+ACDM even the 
worst-case matched IDE model can be distinguished by 
observations. 

We have also shown that it is not always possible 
to construct a physical IDE model which matches the 
growth and expansion histories of our STT models and 
that there is a limit on the strength of deviation from 
GR, beyond which the time derivative of the IDE scalar 
field becomes complex before the present day. 
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Appendix A: Limit on ft in the DGP example 

For the matched IDE/DGP setup studied in Section ILT1 
there is a limit on how close VLi can be set to VLi. The 
differential equation for f2 from [48| which matches the 
DGP and IDE growth histories is, 



b ' 5 ' = 2 S 



1 



— ) fi- (1 

3/3/ 



2Q 2 )n 



where n 2 — 8irG, c 
/? = -(l + 2 )/(l' 
expressed by, 



Q 



is the IDE scalar field, the function 
- f2 2 ) and Q is the coupling function 



3 + 2f 



2^ 

Z H 
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(A2) 



(Al) 



where aH = H. Eq. (jAlj) is quadratic in ft', so to solve 
it for Q we must first solve it for f2', but it is not always 
the case that real roots exist. Using the initial conditions 
specified in [48j the solutions are initially complex for 
tti ~ fli. As f2j is decreased the solutions extend to later 
times but there are no solutions which reach the present 
day for f2j > 0.997. This can be seen in the top left panel 
of Fig. [5] where solutions for values of fii either side of 
this limit are plotted. 



Appendix B: Limit on ft in the scalar-tensor theory 
model and the small uj limit 



For the STT setup of Section Mil the solutions of the 
quadratic Eq. (JTSJ) are not initially complex for f2^ ~ tti 
as they are for Eq. (IA1|) of the DGP setup discussed 
above. A similar solution limit on how close Hi can be 
set to Qi does exist however and depends on uj. In ad- 
dition there is a physical limit which is reached before 
this solution limit and prevents the existence of physi- 
cal IDE counterparts for those STT cases which deviate 
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FIG. 5: Top left: solutions of f2' in the IDE/DGP setup with initial conditions either side of the limit in f2i. Top right: an 
example of ip' 2 becoming negative before the present day in the IDE/STT setup. Bottom left: solutions of Q! in the IDE/STT 
setup continuing beyond the present day, before which ill' 2 has become negative. Bottom right: ip' and Q' limits on f2; as a 
function of ui for the IDE/STT system. 



most greatly from GR. Similar problems have also been 
found in studies of parameterised STT models (l7l. |39|. 

The denominator in Eq. (|16[) can be shown to equal ip' 2 
using Eq's. (jimi3|) . This decreases and reaches zero when 
the universe begins to accelerate and the 4^- term grows 
faster than the f2 term decreases. It can then become 
negative, which would require -0' to be complex and so 
we take this as a physical limit. The top right panel of 



Fig. O shows this happening before the present day for a 
particular choice of oj and O^, while the bottom left panel 
shows that at the same time solutions for fl' still exist. 

We plot both the ip' and CI' limits in the bottom right 
panel of Fig. [5l showing that the smaller the value of uj, 
(and so the greater the deviation from GR), the farther 
fli has to be from Qi . The limit beyond which -0' becomes 
complex shows that it is not possible to find a matched 
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IDE/STT system for u> < 10, contrary to the statement 
in [48j that for any given MG model it is always possible 
to construct a matched IDE model. 

Note that it is possible to finetune ^ to be very small 
and negative so that the t/j' limit is avoided, (too much 
and the STT universe contracts at late-times). Con- 
versely taking <^ small and positive shifts the limit to 
much larger u making it impossible to find a physical 
IDE counterpart for cases with any noticable deviation 
from GR at all. 

The reason that the derivative of the IDE scalar field 
does not become complex for the DGP model can be seen 
from the 'Friedmann' equation of [48| where they define, 



-"0 



(Bl) 



with Cl rc = 0.170. Differentiating this with respect to TV 
and using (l — Cl e~ 3N E~ 2 leads to, 
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This quantity varies from about — 1.50 m at early times 
when E is large, to roughly —2.50 at late times as E — » 1. 
Since f2 m < O at all times we therefore find a condition 
which is true at all times in the IDE/DGP setup, 
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(kc/)') 2 = -3fi m - 2— > 0. 
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